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ABSTRACT 

The r-process in a low temperature environment is explored, in which the neutron emission by 
photodisintegration does not play a role (cold r-process). A semi-analytic neutrino-driven wind model 
is utilized for this purpose. The temperature in a supersonically expanding outflow can quickly drop to 
a few 10 s K, where the (n, 7)- (7, n) equilibrium is never achieved during the heavy r-nuclei synthesis. 
In addition, the neutron capture competes with the /3-decay owing to the low matter density. Despite 
such non-standard physical conditions for the cold r-process, a solar-like r-process abundance curve 
can be reproduced. The cold r-process predicts, however, the low lead production compared to that 
expected in the traditional r-process conditions, which can be a possible explanation for the low lead 
abundances found in a couple of r-process-rich Galactic halo stars. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances — stars: neu- 
tron — supernovae: general 
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1. INTRODUCTION 

It has been generally believed that the high en- 
tropy supernova ejecta (neutrino-driven wind) is the 
likely site for the rapid neutr o n-capture (r- ) pro- 



cess (IWoosley fc Hoffmanl |1992t [Meyer et all Il992t 
IWooslev et alJl994lTakahashi et al.lll994l ). In the previ- 



ous nucleosynthesis calculations related to the neutrino- 
driven wind scenario, the r-process has been examined 
in sufficiently high temperature conditions, such that 
T9 > 1-0 (where Tg is the temperature in units of 10 9 K). 
This physical condition, as well as the high neutron num- 
ber density (N n > 10 20 cm -3 ), has long been consid- 
ered as the physical requireme nts to account for the so- 
lar r-process abun dance curve (|Mathews fc Cowanlll99fl : 
iKratz et al.lll993ft . However, if the neutrino-driven out- 
flow becomes sup ersonic, the temper ature can quickly 
f all below T 9 = 1 (lArcones et al.|[2007h . 

Wan aio et al.l (|2002f ) have examined the r-process in 
such supersonically expanding outflows by introducing 
the "freezeout temperature" Tf (or Tgf in units of 10 9 K) 
that mimics the abrupt wind deceleration by the pre- 
ceding supernova ejecta. In their calculations, the so- 
lar r-process curve around the third abundance peak 
(A = 195) was best reproduced with Tgf « 1.0, and the 
lower Tgf (down to .6) resulted in formin g a lower shifted 
peak. As shown in lArcones et al.1 ((2007) , however, it is 
quite possible that the wind temperature drops to a few 
I0 8 K before deceleration. 

In this Letter, the r-process in such a low temperature 
environment (hereafter "cold r-process") is explored. A 
semi- a nalytic neutrino-d r iven w i nd model (Dunca n et aLl 
19861: lQlan" fc Wooslevl H9M ICardall fc Fulled il997t 



Otsuki et al.H2000tlWanaio 1 et ; aUl200H iThompson et all 
20011 ) is used to obtain the thermodynamic histories of 



outflows (§ 2). Nucleosynthesis calculations in the winds 
are performed with various Tgf including a rather low 
value that has not been considered in previous works 
(§ 3). Fundamental behaviors of the cold r-processing 
are then discussed in some detail. In § 4, a possible abun- 
dance feature in the cold r-process is discussed. Finally, 
implications of this study are presented (§ 5). 



2. NUCLEOSYNTHESIS IN WINDS 

The wind trajectories are obtained using the semi- 
analytic, spherically symmetric, general relativistic 
mode l of neutrino-driven winds (IWanaio et al.1 1200 lL 
2002). The rms average neutrino energies are taken to be 
10, 20, and 30 MeV, for electron, anti-electron, and the 
other flavors of neutrinos, respectively. The mass ejec- 
tion rate at the neutrino sphere M is determined so that 
the wind becomes supersonic through the sonic point. 

The neutron star mass is taken to be 1.4 Mq. The 
radius of neutrino sphere is assumed to be R U {L U ) = 
(i?„o — R v x){L v / L v q) + i?„i as a function of neutrino 
luminosity L„, where i?„o = 15 km, i?„i = 10 km, and 
T„o = 1 x 10 52 ergs s _1 . This mimics the early time evo- 
lution of R y in hydrodynamic simulations (jWooslev et al.l 
1994). The wind trajectories are calculated for 39 con- 
stant L v = (10 — 0.5) x 10 51 erg s _1 with the interval of 
0.25 x 10 51 erg s _1 . In order to link the series of constant 
L v trajectories to realistic time-evolving winds, L v is de- 
fined as L u {t v h) = LyQ^ph/ta) 1 , where i p b is the post 
bounce time and in = 1.0 s. 

The temperature and density in each wind are set to 
be constant when Tg decreases to a given freezeout tem- 
perature Tgf (Fig. 2a), in order to mimic the wind de- 
celeration by the preceding ejecta. In this study, Tgf is 
taken to be 1.3, 1.0, 0.7, 0.4, and 0.1. 

The nucleosynthetic yields in each wind are obtained 
by solving an extensive nuclear reaction network code. 
The network consists of 6300 species between the proton 
and neutron drip lines predicted by a recent mass for- 
mula (HFB-9. lGorielv et al.ll2005l ), all the way from sin- 
gle neutrons a nd protons up to the Z — 110 isotopes (see 
lWanaidl2fjol for more detail). The a-decay chains and 
spontaneous fission processes are taken into account only 
after the freezeout of all other reactions (|Cowan et al.l 
Il999t IWanaio et all |2002| ). Neutron-induced and /?- 
delayed fissions, as well as the contribution of fission frag- 
ments to the lighter nuclei, are neglected. 

Each nucleosynthesis calculation is initiated when the 
temperature decreases to Tg = 9, at which only free nu- 
cleoids exist. The initial compositions are then given by 
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Fig. 1. — Comparison of the mass-averaged yields (line) for Tgf = 
(a) 1.3, (b) 1 0, (c) 0.7, (d) 0.4, an d (e) 0.1 with the solar r-process 
abundances (Kappclcr et al. 1989, dots) scaled at the height of the 
third r-process peak (A «s 195), as functions of mass number. Color 
and black lines are the yields before and after a-decay, respectively. 

the initial electron fraction Y e (number of proton per 
nucleon). In order to mimic the hydrod ynamic results 
(|Wooslev et al.lll994HArcones et al.ll2007t ). Y e is assumed 
to be Y e (L v ) = (Y e o - Y el ){L u /L u0 ) + Y el . Y e0 and Y eX 
are taken to be 0.50 and 0.15, respectively. The latter is 
rather s maller than that in hydrodynamic results (e.g., 
~ 0.35 in Woosley et al. 1994). This extreme assumption 
is to obtain the heavy r-nuclei abundances with a typical 
neutron star mass (1.4 Mq) for the current purpose. 

3. HOT r-PROCESS VS. COLD r-PROCESS 

The nucleosynthetic yields in each Tgf case are mass- 
averaged over the 39 wind trajectories weighted by 
M(L u )At v \,. Figure 1 compares the mass-averaged yields 
with the solar r-process abundances as functions of 




t (s) 



(b) N n , (c) T„ 7 /r 7n , (d) 
s -1 as functions of time. 



Fig. 2. — Time variations of (a) Tg 
Tp/Tn-y, and (e) 5 a for L„ = 2 X 10 51 ei 
The assigned color in each Tg{ case is as the same as in Figure 1. 
The times corresponding to the second (A = 130) and third (A = 
195) peak formation, n-exhaustion, and frcczeout are marked by 
open and filled triangles, and open and filled circles, respectively. 
The dotted lines in the bottom panel show 5^ (see text). 

atomic mass number. We find that, for Tgf = 1.3 
(Fig. la), the third-peak abundances shift to the high 
side from A = 195. In addition, the abundance curve is 
unacceptably jagged compared to the solar r-process pat- 
tern. For Tgf = 1.0 (Fig. lb), the third-peak abundances 
fit the solar r-process pattern quite well. For Tgf = 0.7 
(Fig. lc), however, the third-peak abundances slightly 
shift to the low side from A — 195. 

This can be understood within the classical picture of 
r-process, in which the (n, "/)-("/, n) equilibrium is justi- 
fied. That is, for lower Tgf, the r-process path locates 
at the more neutron-rich side on the nuclide chart. As a 
result, the freezeout at N = 126 takes place at lower Z 
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Fig. 3. — Comparison of the mass-averaged yields (colored 
lines) with abundances for CS 31082-001 ( circles, Iffill et al.ll2002l ; 
llvarsson, S. et al.l [20031: IPlez et al.ll200H) and the solar r-proccss 
curve (black lines. iBurris et all 120001 ) scaled at the Eu (Z = 63) 
value, as functions of atomic number. The assigned color in each 
Tgf case is as the same as in Figure 1 . The vertical scale for the up- 
permost set is true, and the others are scaled downward for display 
purposes. 

(i.e., lower A). Among the above three cases, Tgf = 1.0 
can be regarded as a reasonable choice to reproduce the 
heavy r-process nuclei in the solar system. We find, how- 
ever, that this interpretation does not hold in the lower 
Tgf cases. The Tgf = 0.4 (Fig. Id) case reproduces the 
third-peak abundances even better than the Tgf = 0.7 
case. Moreover, for Tgf = 0.1 (Fig. le), the third-peak 
abundances nicely fit the solar r-process curve. This is a 
consequence of the cold r-processing as discussed below. 

Figure 2 shows some key physical quantities as func- 
tions of time (t = at Tg = 9) in the L v = 2 x 
10 51 erg s" 1 case. The asymptotic entropy per nucleon 
in this wind is only 92 (in units of the Boltzmann con- 
stant). However, the (artificial) low initial Y e (= 0.22) 
leads to a high neutron-to-seed ratio Y n /Yh (= 109) at 
the beginning of the r-process (Tg = 2.5). Here, Y n and 
Yh are the abundances of free neutrons and the heavy nu- 
clei [Z > 2), respectively. At Tg = 2.5, the abundance- 
averaged mass number of the heavy nuclei is 103. This is 
thus taken to be representative of the third-peak forming 
winds. On each line in Figure 2, open and filled triangles 
denote the times when Y n /Yh decreases to 80 and 10, 
respectively. The former and latter approximately corre- 
spond to the epochs at which the second (A = 130) and 
third (A = 195) peaks form, respectively. Open circles 
denote the time when Y n /Yh decreases to unity, which 
is referred to as "n-exhaustion" . No heavier nuclei are 
synthesized beyond this stage, but the rearrangement of 
the local abundance distribution continues. 

The ratios r„ 7 /r 7n and ra/T„ 7 in each case are shown 
in Figures 2c and 2d, where t„ 7 , r 7 „, and are the 
abundance-averaged mean lifetimes of (71,7), (7,rt), and 
/3-dec ay for Z > 2 nuclei , respectively (see eqs. (3) and 
(4) in Wa naio et al.|[2004T ). These represent the lifetimes 
of the dominant species at a given time. The condition 
Tp/Tn-y = 1 is referred to as "freezeout" (filled circles] 
Fig. 2). The abundance curve fixes at this time. Note 
that the classical r-process is characterized by the con- 



ditions t„ 7 /t 7 „ « 1 and r^/r„ ^ 1. 

As can be seen in Figure 2c, the second peak forms 
at t = 0.08 — 0.1 s in all cases, where t„ 7 /t 7 „ ks 1. 
This means that the second peak forms under the con- 
dition of the (rt, 7)-(7, n) equilibrium. For Tgf = 1.3, 
this condition holds during the whole stage of the re- 
process. For Tgf = 1.0, t„ 7 /t 7 „ slightly decreases after 
rt-exhaustion and a quasi (n,7)-(7,n) equilibrium con- 
tinues until freezeout. For Tgf = 0.7, t„ 7 /t 7 „ slightly 
decreases after the second-peak formation and a quasi 
equilibrium continues by the third-peak formation. In 
these cases, therefore, the temperature is enough high 
to maintain the (quasi) (n,7)-(7,n) equilibrium during 
the major r-process phase. Hereafter, the Tgf = 1.0 case 
is taken to be representative of the r-process in a high 
temperature environment (hereafter "hot r-process"). 

In contrast, in lower Tgf cases, t„ 7 /t 7 „ quickly drops 
below unity and the (n, 7)-(7, rt) equilibrium is never 
achieved after the second-peak formation. This is due 
to the quickly decreasing Tg and N n , in which the N n - 
Tg condition a bruptly falls off the "waiting point valid- 
ity b o undary" (I Cameron et al. [l983: Mathew s fc Cowanl 



119901 IGorielv fc Arnouldl 119961 ). The temperature does 
not play any roles in such an environment, which is thus 
designated as the cold r-process. Hereafter, the Tgf = 0.1 
case is taken to be representative of the cold r-process. 
One may consider that, without photodisintegration, the 
r-process path approaches the neutron-drip line and the 
third peak shifts to the rather low side from A = 195. 

The reason for this misunderstanding can be found in 
Figure 2d, where the lifetimes of /3-decay and of neutron 
capture are compared. The hot r-process satisfies the 
classical r-process condition 173 /t„ ^ 1. In the cold r- 
process, however, 773 is only a few times larger than t„, 
in which the neutron capture competes with the /3-decay. 
This is due to the quickly decreasing N n in the supersoni- 
cally expanding wind. Consequently, the nucleosynthetic 
path of the cold r-process is pushed back and locates at 
a similar position to that of the hot r-process. 

This is clearly seen in Figure 2e that shows the re- 
process path in each case in terms of the neutron separa- 
tion energy. Here, Sinj'l (two-neutron separation energy 
divided by two) is abundance-averaged for the Z > 2 
nuclei (S a ) , which approximately represents the nucle- 
osynthetic path at a given time 1 . That is, S a = is the 
neutron-drip line, while a higher S a locates at closer to 
the /3-stability. The path obtained fr om the (n, 7M7, n) 
equil ibrium condition (S®, eq. (3) in IGorielv fc Arnouldl 
1996) for each is also plotted, which poorly predicts the 
real path for the cold r-process. The cold r-process takes 
a significantly higher S a than S% owing to a push-back by 
/3-decay. As a result, both the cold and hot r-processes 
trace similar S a histories. Notably, the freezeout in the 
Tgf = 1.0 and 0.1 cases takes place at similar S a values. 
As a consequence, similar r-process curves appear as can 
be seen in Figures lb and le. Note that, for Tgf = 0.1, 
freezeout occurs at Tg « 0.2. Hence, any Tgf < 0.2 cases 
would result in the identical abundance curve. 



4. ABUNDANCE FEATURE IN THE COLD r-PROCESS 

1 The animations for Tgf = 1.0 and 0.1 are available from 
http/ /supernova. astron.s. u-tokyo.ac.jp/~wanajo/research. html. 
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TABLE 1 
Mass Fraction (in units of 10 -2 ) 



T 9f 


Zr 


Ba 


Eu 


Os 


Ir 


195 pt 


Pb 


Th+U 


1.3 


7.58 


0.56 


0.14 


0.19 


2.59 


1.86 


1.67 


1.40 


1.0 


7.97 


1.89 


0.20 


2.04 


3.24 


1.69 


2.06 


1.27 


0.7 


8.04 


1.49 


0.31 


2.90 


3.51 


1.05 


1.45 


1.26 


0.4 


8.02 


1.32 


0.32 


2.71 


3.49 


1.25 


0.97 


1.31 


0.1 


7.97 


1.32 


0.36 


2.34 


2.97 


1.58 


0.92 


1.25 



The abundance differences between the hot and cold r- 
processes become further small when plotted as functions 
of atomic number. Figure 3 compares the mass-averaged 
abundances with the solar r-process curve and the mea- 
sured elements in the r-process-rich Galactic halo star 
CS 31082-001. All cases, except for T 9f = 1.3, are in 
reasonable agreement with the solar r-process and mea- 
sured patterns. This may support the robustness of the 
abundance curves found in the r-process-rich stars. Note 
that the total mass of the r-process nuclei (Z > 90) is 
the same for all cases (« 2.9 x 10~ 4 M©), since the seed 
abundances (A ~ 100) are already produced at Tg ~ 2.5 
(> Tgf). The mass fractions of some measurable elements 
(and l95 Pt) are presented in Table 1, which also shows a 
small variation between the hot and cold r-processes (up 
to a factor of two, except for Tgf = 1.3). These differences 
are small enough, when considering the uncertainties in 
the nuclear data as well as in the outflow conditions. 

We can see, however, a distinctive feature in the abun- 
dance curves around A = 200—230 (before a-decay; black 
lines in Fig. 1) between the hot and cold r-processes. 
This is due to their different freezeout histories (Fig. 2e) . 
The hot r-process (red line) follows a low S a (« 2 MeV) 
path at n-exhaustion. At this S a path, a trough at 
A w 200 appears in the abundance curve (before ot- 
decay, Fig. lb). This is a consequence that the nuclei 
at N w 130 — 140 are predicted to be unstable against 
neutron capture (s ee the iso-S^n/^ curves in Fig. 6, 
IWanaio et alj [2004 ). On the other hand, the cold r- 
process (blue line) takes a high S a (~ 3 MeV) path 
at n-exhaustion (Fig. 2e). Therefore, the trough shifts 
to A « 210. The trough is further eroded by neutron 
capture without push-back by photodisintegration. As 
a consequence, the parent nuclei of Pb (A = 210 — 231 
and 234; Fig. le) result in being significantly deficient. 
This may lead to a visible difference in the Pb production 
(Fig. 3 and Table 1), despite the similar Th+U values 2 . 

5. IMPLICATIONS 

In this Letter, the r-process in a low temperature en- 
vironment (Tgf ^0.1; cold r-process) was investigated, 
which is characterized by the conditions t„ 7 /t 7 „ <C 1 

2 There is a local problem in the odd-even behaviors for A > 240 
(Fig. 1) when the current nuclear data set is utilized. This leads to 
an unusual pattern among 232 Th, 235 U, and 238 U, which however 



and Tg/r n7 > 1. Despite the fundamental difference of 
the cold r-process compared to the traditional (hot) re- 
process, a solar-like r-process curve can be reproduced 
quite well. These physical conditions are realized in 
the supersonically expanding neutrino-driven outflows. 
In particular, the winds from a low-mass progenitor 
that has a steep core density gradient easily expand 
to a large distance (and thus to low T gf) before col- 
liding to the preceding supernova ejecta (jArcones et al.1 
l2007f) . Therefore, the low- mass end of supernova pro- 
genitors (e.g., 8 — 12 M©) may be the most likely site 
of the cold r-process. The low-mass supernovae have 
been also suggested to be the r-process site from Galac- 
tic chemical evolution stu dies ([Mathews fc Cowanlll990l : 
llshimaru fc Wana id [19991 ). 

Among the six r-process-rich Galactic halo star s cur- 
rently reported, CS 31082 -001 dPlez et alj I2004D and 
HE 1523-0901 (upper limit. iFrebel et al.ll2007f) appear to 
have the low Pb abundances compared to the scaled solar 
r-process curve. Th e upper limit of Pb for BD +17° 3248 
(|Cowan et al.l [2002;') is still consistent to the scaled so- 
lar va lue. The othe r three, CS 22892- 052 ([Sneden et all 
[20031). HP 115444 dWestin et al.l[200l . and HD 221170 
(|Ivans et al.1 12006[ ) have the measured Pb abundances 
consistent to the (high) scaled solar value. The varia- 
tion of the Pb production may be attributed to the dif- 
ferent types of (cold and hot) r-processes. The low Pb 
abundance in CS 31082-001 has been in fact a most wor- 
risome aspe ct, when the U/ Th value is applied for the 
age dating (|Plez et al.ll2004| ). This is hardly explained 
from the classical (or hot) r-process view, in which a low 
S a (« 2 MeV) event is inevitable for the third-peak for- 
mation. In this regard, the cold r-process might enables 
us to consistently treat the Pb, Th, and U abundances 
for the cosmochronology. 

Finally, it is noted that the cold r-process brings us an 
additional challenge to the theoretical r-process study. In 
most of previous r-process calculations, the uncertainties 
in neutron-capture rates have not been considered seri- 
ously. This is reasonable for the hot r-process, in which 
uncertainties even by a factor of ten would not cause a 
serious problem (see Fig. 2d). In the cold r-process, how- 
ever, the neutron capture competes with the /3-decay, and 
thus accurate neutron-capture (as well as /3-decay) rates 
will be required to predict reliable r-process abundances. 
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does not affect the initial Pb abundance. Therefore, only the sum 
of Th and U are presented. 
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